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In usual particle models, sterile neutrinos can account for the dark matter of the Universe only 
if they have masses in the keV range and are warm dark matter. Stringent cosmological and 
astrophysical bounds, in particular imposed by X-ray observations, apply to them. We point out 
that in a particular variation of the inert doublet model, sterile neutrinos can account for the dark 
matter in the Universe and may be either cold or warm dark matter candidates, even for masses 
much larger than the keV range. These Inert-Sterile neutrinos, produced non-thermally in the early 
^\ • Universe, would be stable and have very small couplings to Standard Model particles, rendering 

very difficult their detection in either direct or indirect dark matter searches. They could be, in 
, principle, revealed in colliders by discovering other particles in the model. 
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. Since the first indications of the existence of dark matter more than seven decades ^o_ll|, many different strong 
£NJ ' pieces of evidence in its favor have been accumulated (for reviews, see eg. Refs. @, U Hf). The presence of dark 
C\| , matter has been revealed so far through its gravitational effects. Much effort is being devoted to the detection of dark 
matter annihilation or decay products or the scattering of dark matter particles off nuclei. However dark matter may 
consist of particles which will not be revealed (at least in the near future) in this type of searches. We provide here an 
O .■ example of a dark matter candidate found in a simple extension of the Standard Model (SM), whose nature could be 
indirectly proven only through the discovery and study in colliders of other non-standard particles predicted within 
the model. The dark matter particle candidate we study here is a sterile neutrino with mass in the tens of keV to the 
tens of GeV range and produced non-thermally in the early Universe, which can be either warm dark matter (WDM) 
!— i, or cold dark matter (CDM). 

One or more gauge singlet right-handed (sterile) neutrinos are included in simple extensions of the SM which can 
easily accommodate neutrino oscillation data [H, & & 9 • These data show that at least two of the active neutrinos 
have a non-zero mass. In many models sterile neutrinos are the right-handed Dirac mass partners of the active 
f-^. | neutrinos. In some see-saw-inspired models, sterile neutrinos have large Majorana masses, which leads to three light 
(mostly active) neutrinos and several heavier (mostly sterile) neutrinos, the lightest of which is an attractive dark 
matter candidate [8[. Since this candidate necessarily decays into a light neutrino and a photon, to constitute the 
^vq dark matter its lifetime must be much longer that the age of the Universe. Thus, this dark matter candidate might be 
y—i ' detected through the photons produced in its decay in the dark halos of galaxies. Moreover, to have the required dark 
0^ . matter relic density, the lightest sterile neutrino must usually have a mass in the keV range, although this depends 
■ on the mechanism through which sterile neutrinos are produced in the early Universe. 

Relic sterile neutrinos with only standard model interactions are produced in the early Universe through active- 
sterile neutrino oscillations. Sterile neutrinos produced through non-resonant oscillations @, 0, [1] must have masses 
M s in the keV range to account for the whole of the dark matter and are WDM. Through a combination of X-ray 
and structure formation constraints, an upper bound M s < 3 — 4 keV has been obtained 

@, HUM EH (see, 

however, Ref. [l2j for a very recent weak hint of a possible signal). Lyman-a forest data has been used to impose 
the lower bound of M s > 5.6 keV [l3[ (see also Refs. [14j, [15J for previous bounds) on non- resonantly produced sterile 
neutrinos, or the revised limit of M s > 8 keV obtained by a new analysis which combined with the previous upper 
bounds would exclude non-resonantly produced dark matter sterile neutrinos. Even disregarding the controversial 
Lyman-a bounds, the mass range allowed for these neutrinos is very restricted because there is an independent lower 
bound M s > 1.8 keV (T3, HH derived from the analysis of phase space density evolution of dwarf spheroidal galaxies. 
In general, these bounds do not consider the possibility of a very large lepton asymmetry. In the presence of a 
large lepton asymmetry C = (n Ue — n&^/s > 10~ 6 , where n„ c and np e are the number densities of neutrinos and 
antineutrinos and s is the entropy density in the Universe, sterile neutrinos may be produced in the early Universe 
through resonant oscillations (19 . l20j . Considering the upper limit of the lepton asymmetry imposed by Big Bang 
Nucleosynthesis (BBN) £ < 2.5 x 10~ 3 (2(J, the range 1 keV < M s < 50 keV is in principle allowed for sterile 
neutrino dark matter [l7l l20l. l2l| . In slightly more complicated models, sterile neutrino dark matter may be produced 
as decay products of, for example, a heavy singlet scalar [22, 23], or may not completely thermalize as in low reheating 
temperature scenarios (24|. Yet, in all these models the X-ray constraints are important. 

Here, we consider a small variation of the SM in which the lightest sterile neutrino is stable (hence it does not 
produce photons as decay products) and may constitute all of the dark matter. We study a variation of the Inert 
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Doublet Model [25|, [26| (in itself an extension of the model in Ref. [22]]). In this model one scalar doublet, r\ = (rj + , rjo) 
and three sterile neutrinos, which we call Inert-Sterile neutrinos, iVj with i = 1,2, 3, odd under a new parity Z 2 , are 
added to the SM. All the particles in the SM are even under the additional Z 2 symmetry. These assignments make 
the new particles "inert" because their couplings to the SM particles are very limited. The leptonic Yukawa couplings 
in this model are 



+ </>\)K + h i3 (virf - IjV^Nj + h.c. 



(1) 



Here <f> = (0 + ,0°) is the SM scalar doublet field, and L = (i^,Zj) are the SM lepton fields. Under the extended 
electroweak symmetry SU(2)l x U(1)y x Z 2 , the fields 77, N, </> and L are in the (2,1/2;-), (1,0;-), (2,1/2;+) and 
(2,-1/2;+) representations respectively. The inert and the standard doublet scalar also couple through the scalar 
potential [26j], 



V = nl&Q 



A x ($t$)2 + \ 2 (rf v ) 2 + A 3 ($ t $)(7 ? t r;) + A 4 ($ t ?7)(?7 t *) 
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In particular the last quartic coupling provides the mass splitting between the two physical inert neutral scalar particles 
t]h = v2 Im(?7 ) and tjl = V% Re(?7 ) |2a,|2a|, which are the heaviest and the lightest for A5 < (otherwise the two 
would be exchanged) 
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The masses of the inert scalar bosons are, 
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= fi 2 2 + X 3 v 2 /2 

= fi 2 2 + (A 3 + A 4 - X 5 )v 2 /2 = n% + (A L - 2X 5 )v 2 /2 
= A + (A3 + A 4 + X 5 )v 2 /2 = nj + X L v 2 /2 . 



(4) 



Here v/y/2 = 174 GeV is the vacuum expectation value (VEV) of the SM Higgs field (the inert scalar does not acquire 
a VEV), m v ± is the mass of the charged scalars, A5 has been chosen to be real and we define A^ = A 3 + A4 + A 5 . 
The only other terms in the Lagrangian allowed by the Z2 symmetry are Majorana mass terms for the Inert-Sterile 
neutrinos, 
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~-M i N i N l + h.c. 



(5) 



The Zi symmetry forbids sterile-active neutrino mixings. The N^s are not the Dirac mass partner of the Vi as in 
usual extensions of the SM and active neutrino Majorana masses are generated at one-loop level. The active neutrino 
mass matrix elements are 25] 



{M v )ij = ^2 h ik h 3k 
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We will assume in what follows that m ri 



is of the order of 100 GeV and m nL of the order of tens of GeV, thus the 



first term in Eq. [6] is dominant. 

The Z2 parity im plie s that the lightest inert particle is stable and thus a good dark matter candidate. Both the 
lightest inert scalar |25l . l2rl [H, [29|, |3(| HH, [H, [33| and the lightest sterile neutrino [HI, [H, H|[ could be dark matter 
candidates. We will assume the second possibility. 

In Refs. [34I. l3a] it was assumed that the mass difference between tjl and t]h is small, i.e. the coupling A5 is very 
small. In this case, in order to generate the observed active neutrino masses, the hij couplings cannot be very small. 
In addition, it was assumed that mo = {m 2 + m 2 )/2 > Mi,M 2 ,M 3 and the lightest Ni is produced thermally. 



Under these assumptions, Ref. [34| found that the lightest Inert-Sterile neutrino can be CDM and account for the 
whole of the dark matter if its mass is in the range 7 GeV to 300 GeV. 

Here we will explore a range of values of the coupling constants different from those previously considered, namely A5 
not very small and Yukawa couplings small enough to ensure that the sterile neutrinos Ni are never in equilibrium 
in the early Universe. We will not study the flavor structure of the couplings hij, but only their order of magnitude. 
We call generically hx, /12, hi the couplings of Ni, N2 and A/3, respectively. We assume a hierarchy in the couplings, 
with h% < h,2 — h^. We also assume that only the lightest sterile neutrino, which we take to be N%, is lighter than the 
lightest inert scalar tjl and hence, it is the dark matter candidate. The tjl particles are produced thermally in the 
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FIG. 1: Diagram for fj, — » e7 transition in the variation of the inert doublet model we consider here. 

early Universe and decouple when they are non-relativistic. The subsequent late decay of the t]l produces the Inert- 
Sterile Ni relic particles that now constitute the dark matter. In this scenario, depending on the mass, abundance 
and lifetime of T)l, the N% can be either CDM or WDM and account for the whole of the dark matter with mass in 
the range ~few keV to tens of GeV. We will show that all requirements on the model are fulfilled: active neutrino 
masses of the right order of magnitude are obtained, the upper bound on the hij from fi — > ej is easy to fulfill, all Ni 
producing reactions in the early Universe are out of equilibrium and the necessary relic density and decay rate of t]l 
for different values of the t]l and N\ masses are obtained, while respecting all the collider and other bounds imposed 
on the model. 

Let us see first how large the Yukawa couplings h must be to get reasonable values for the active neutrino masses, 



10 1 eV. Using Eq. [SI and assuming that tjl is significantly lighter than t]h, that M 2l 3 is of the same 



order of magnitude but larger than m VH and that the contributions of iV2 and N3 are dominant, we get 
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Notice that when m TIH is large with respect to m rjL , Eq. implies that m rlH 
larger than but similar to m riH , the logarithm in Eq. [T] is close to 1, thus 
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Lepton flavor transitions like the fi — > e'y process in Fig. [T] occur in this model. The branching ratio, B^ c 



^u— te^l^iM— >euv in the inert doublet model is [34l. |3(: 
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where a is the fine structure constant and Gf is the Fermi constant. For M\ 
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experimental upper bound on the branching ratio, B(fi 

h 2 , 3 < 2 x 10 



1/12, whereas for Mi < m n ± it is ^(0) 
> ei) < 1.2 x 10" 11 [37], implies 
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for hi « ft.2,3- 

Let us now see how small the couplings hij must be in order for the Ni to never be in equilibrium in the early 
Universe. The upper bounds are particularly important for N2 and N3, whose generic couplings, I12 and /13, are 
larger than the coupling h\ of N\. The Ni can be produced through the reactions in Fig. [2j i.e. two to two reactions 
LL —>■ NiNi mediated by any of the four physical inert particles r]H,r]L,r] + ,rj~ , which here we call now generically 77, 



or 7777 — * NiNi mediated by L. The Ni could in principle be produced through the decay 77 
decay of rjL — * JVj. The production rate for N 2 , for example, is 



NiL and the inverse 
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FIG. 2: Production processes of Inert-Sterile neutrinos in the early Universe. 

where riff is the N2 equilibrium number density which appears in the equation as a normalization factor, and ul and 
n„ are the number densities of the standard leptons and the inert scalars respectively, at the temperature considered. 

Eq.[Tl]is derived from the Boltzmann equation for the production of Ni (i — 1, 2, 3) in the process ab — ► A^c, where 
a, b and c are particles and we assume that only the initial particles, a and 6, have initially a non zero particle density. 
If the particles a and b are in equilibrium, assuming Maxwell-Botzmann density distributions, the time evolution of 
the number density depends on the number densities of particles a and b in the following manner (see e.g. Eqs. 
(5.8) and (5.23) of Chap. 5 of Ref. [H) 

— — + 3Hn Ni = n a n b < a ab ^ NiC \v\ > . (12) 
As usual, it is convenient to change variables to Y = tinJ s and x = m^jT (see Eq. (5.16) of Ref. [38|) and obtain 



V* n a n b < <7 ab _,iv iC |u| > (13) 

a,b,c 



dx H X S 



Now, dividing and multiplying the right-hand side of Eq. [T^l by normalizing function, one gets 

_x_dY_ = 

Y e i dx H v ' 

where is defined as in Eq. QT] above. Eq. [14] is equivalent in this case to Eq. (5.26) of Ref. (38| and shows that 
Yat is never significantly different from zero if Y^ i /H < 1. 

Following Refs. [35^ |4(| and making use of Refs. [4l|, |42j, for relativistic 77, Ni and L the thermal averaged cross 
section of LL — ► NN and of r\r\ — > A^A^ are approximately given by 

< <™ >£T-^jv« - 0-7 x 10- 1 1| , (15) 

h 4 ( AT 2 \ 

< av > m ^ NiNi =s 3 x 10 ^ In ^ M a + 7n 2 J ' ( 16 ) 
which show that the process r\r\ — > A/jA^ is dominant and 

4T : 

The production is out of equilibrium if the rate is smaller than the expansion rate of the Universe, H 



F Ni ~ IV^^ ~ 0.7 x 10- 1 hf T In ( M ~ ~ m2 ) . (17) 



r Ni <H = im y fg:^— , (18) 

where is the number of degrees of freedom and Mp\ is the Planck mass. Since the right-hand side of Eq. [TS] 
decreases faster than the left-hand side for decreasing T, if the condition is fulfilled for the smallest T value in the 
range considered, i.e. the smallest T for which all the particles involved in the production are relativistic, then it is 
fulfilled for all larger T. 

At high temperatures T > M^,a — rn VH we need to write the condition in Eq. [TH]at T ~ Mk ~ m VH . Thus, the 
production of relativistic ^2,3 is out of equilibrium at T > M-i^ ~ m riH if 

1/4 
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Since we are assuming Mi < m VL << M2.3, the condition in Eq. [18] for relativistic N% and r\h must be taken at 
T ~ , thus the production of relativistic Ni from relativistic 77^ is out of equilibrium if 

hi < 2 x 10- 4 f -^-V" ( ^^V /4 . (20) 

1 V 106.75/ VlOGeW v ' 

At temperatures in the range M2.3 > T > m^, in which the ^2,3 are non-relativistic (but 77^ and L are relativistic), 
the relevant thermal average cross sections for A^ production are approximately 

, hf ( 2MA , « 

< ov > LL ^ NiNi c~ 0.8 x 10- 2 ^ cxp {--f 1 ) > (21) 

1 ( 2AL, 



< av > nn ^ NtNt 2 x lO^ 1 ^ exp I ^ ] . (22) 



The production is again dominated by the r\Lf]L — ► NiNi process, thus 



T>5/2 / »r 

-1 1,4 1 / 



IV, ~ r ww ^^ 4 ~ 0.7 x HT 1 /if — ^ exp (-jr) ■ (23) 

Because this rate decreases faster than H , if < 77 is fulfilled at T — A/2.3 where T is maximum within the T 
interval, the process will be out of equilibrium for lower values of T, thus we obtain 



h 2 , 3 < 3 x 10" 



■*' Y /8 (J^) 1/4 (24) 
J llOOGevj ■ [ ] 



106.75/ VlOOGeV; 



For still lower temperatures T < m v ±,o, for which all the inert bosons are non relativistic but the Ni are relativistic, 
we need to verify that the N\ are not produced thermally (recall we are assuming that m^'° > Mi). In this case 

3 4 M? 

< av > vv ^ NlNl ~ — h x —£ , (25) 

and 

~ T, p ^ NlNl = 10- 2 hf Ml exp (-^A . (26) 

This rate decreases faster than H as T decreases, thus if Tn 1 / H < 1 at the highest temperature in the range considered, 
T = m v , the condition is fulfilled at any lower T. Thus, 

hl < 3 x 10- 2 f^) 1/S (^) 3/4 fM£XV /2 . ( 27 ) 
V 106.75/ VlOGeV/ V Afi / V 7 

After considering all the required upper bounds on the hij Yukawa couplings, we conclude that Eq. fTS] provides 
the most restrictive upper bound on the Yukawa couplings of the heaviest inert sterile neutrinos, /i2,3, and they are 
compatible with the value assigned to h 2t 3 in Eq. [7] which is necessary to account for the active neutrino masses. The 
most restrictive bound on hi, the Yukawa couplings of the lightest Inert-Sterile neutrino TVj., will be given in Eq. [29] 
below and is derived from our requirement of a long enough lifetime of the lightest inert bosons rjL into Ni. 

Let us now consider the decays of the 77 ± and t]h into Inert-Sterile neutrinos. If > to° + raw, then the 
process 77 ± — > 77° + W can occur. The branching ratio of the decay mode rf 1 — ► NiL with respect to the dominant 
r / ± ~~ * V° + W mode is proportional to the ratio of the couplings hf/g w , where gw is the weak coupling. Using 
the value of /i2,3 necessary to produce the active neutrino masses, given in Eq. [H] with |As| ~ 0.2 for example, 
hf/g w ~ 10~ 10 (M^/lOO GeV), which is negligible. Thus, the heavier Inert-Sterile neutrinos N 2} 3 are not produced 
in the decays of the inert charged bosons. Neither the lightest Inert-Sterile neutrino is produced in these decays, 
since hi << ft-2,3- If, instead m^± < m„o + raw , the 3-body decay 77^ — > 77 + L + L dominates the decay of 77^; the 
branching ratio of 77 * — > NiL then goes as hf/g w ~ 10" 11 (Mj/10QGeV) for the heavier Inert-Sterile neutrinos. The 
branching ratio is even smaller for TVi. Again, the decay of the charged inert bosons into the Inert-Sterile neutrinos 
Ni is negligible. For the decays of the heavier neutral inert boson 77^ the same arguments apply but changing the 
W's by Z's. Thus the Inert-Sterile neutrinos are not produced in the decays of r/ ± and 77^. 



FIG. 3: Dominant r/L annihilation channels into Standard Model fermions /. 



We need to insure that 77^, the lightest inert scalar particle, is produced thermally in the early Universe and that 
it is in equilibrium before decoupling while it is already non-relativistic, at freeze-out, T/. . < m VL . The dominant 
processes that maintain the 77^ particles in equilibrium depend on the couplings of t]l with the SM particles. The rj^ 
gauge couplings and its couplings in the scalar potential are the same that occur in the inert doublet model in the 
absence of sterile neutrinos. Using the same couplings, in Ref. [26|, [2^,[32| T)l with mass in the GeV range are found to 
be good dark matter candidates. We want instead that the tjl decay into the lightest inert sterile neutrino N±, which 
constitutes the dark matter now. After the 771, particles decay through the process t\l - > NiVi, there is one N\ per 
each rjL. In order for N\ to account for the whole of the dark matter, the number density of t]l at their de coup ling 
must be larger for the case considered here than in the scenarios in which they constitute the dark matter [2(| [U H2] • 
The number density that is needed for non-relativistic N± to be the dark matter at present, must be the same 
relic number density n VL the 77^ should have at present had they not decayed. Thus the relic density of N\ is now 
n^Mi = n VL Mi and 



n Nl h 2 =VL r]L h 2 ( — i ) , (28) 



Mi 

where fi^/i 2 is the relic density the 771, would have at present if they were stable. When the N\ can be either CDM or 
WDM we require the N\ density to be that of the observed relic density of dark matter f2 DM ft, 2 = 0.1099 ±0.0062 43]. 
If the Ni are instead hot dark matter (HDM) we should impose the upper bound flN 1 h 2 < 0.014 = ^HDM-max^ 2 (the 
95% CL on the relic density of light neutrinos) [43| . 

If m nL > raw-, the 77^ annihilate efficiently into two W bosons and their relic density is too small even to constitute 
the bulk of the dark matter, thus we are not interested in this mass range. When m r]L < mw, the processes in 
Fig. [3] and their inverse processes keep ?7l in equilibrium. The lightest scalar T)l coannihilates with the heaviest 
inert scalar partner tjh- The coannihilation r\nr\h —* ff into SM fermions / is mediated by the Z boson and its 
cross section depends on the mass spittling A = m VH — m VL which in turn, depends on A5 (see Eq. [3]). The tjl also 
coannihilates with 77^, via W ± exchange, with a cross section which depends on the mass split between them. The 
process r\hT\L — *■ // via Higgs exchange also keeps 77^ in equilibrium, and in the particular range of masses we explore 
below is the dominant process. We use the public code MicrOMEGAs [44( to compute the 77^ relic density. 

The decay 77^ — > N\L must happen after the 77^ freeze-out at Tf = m VL /xf where Xf is in the 20 to 30 range. 
Thus, the decay rate must be T 71 ^,n 1 l — h\m riL /\Qi: < H for T > Td eC ay and T 11 ^n 1 l ~ H for T = Tdccay with 
Tdecay < Tf, Q .. These conditions lead to the most stringent bound on hi 



\xf \ 10 GeV/ V 10.75/ V ' 



20 \ / m VL y/a / g* \V 4 
x f J \ 10 GeV/ V 10.75. 

Note that this bound on h\ is consistent with the previous requirements. 

We can now show that the Inert-Sterile neutrinos produced in this model may be either WDM or CDM, which are 
characterized by the free-streaming length A / s [45|, 14a ] 




Af s = 2r/jEQ (1 + ^eq) In W 1 H — jTi — \ v? ^ T\ — I v? • ( 30 ) 



Here the subscript EQ denotes matter-radiation equality and r = a(t)p(t)/Mi, where a(t) and p(t) are the scale 
factor of the Universe and the dark matter particle characteristic momentum at time t respectively. As the Universe 
expands, the ratio r remains constant. At the time of matter-radiation equality, A/ s must be 0.1 Mpc [47j for WDM, 
which fixes r ~ 10~ 7 . At the moment of decay of the 77^ (we make the approximation of instantaneous decays) the 
scale factor of the Universe is a ~ T /Td eC ay, where T a is the photon temperature today, and the momentum of the 
relativistic N\ decay products is m VL /2. Thus, r ~ T m VL / (2 Td GC ay Mi)- Therefore, r = 10 -7 fixes the mass of Ni 
to be 
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FIG. 4: In both panels, Mm gg s = 160 GeV, m, B = 125 GeV and rn r] ± — 130 GeV. Upper panel: The shaded areas correspond to 
forbidden values of Xl from vacuum stability (cross-hatched violet region) and perturbativity (shaded gray region) arguments. 
Lower panel: From top to bottom, the blue, red and green colored narrow strips show the regions where Ni would have the 
right dark matter density for the corresponding values of Xl- The unshaded background region corresponds to the range in 
Eq. 1321 where the Ni may constitute WDM (above it, Ni can only be CDM and below it, only HDM). For any particular 
value of Td cca y between 5 MeV (upper boundary of unshaded region) and m VL /xf (lower boundary of unshaded region, which 
depends on Xl through Xf) there is one value of Mi, given by Eg. 1311 within the unshaded background region for which Ni 
would be WDM (and it would be CDM for all larger values of Mi and HDM for all smaller ones). In order for Ni to be allowed 
as HDM, its mass must be at least a factor of f2DM/i 2 /^HDM-max/i 2 = 0.1099/0.014 ~ 8 smaller than that corresponding to 
the center of the colored bands for a given m VL . 



Given a particular Td eC ay, Eq. [3T] provides the Ni mass for which the Ni would constitute WDM. Heavier N\ (smaller 
Af s ) would be CDM and lighter ones (larger Af s ) HDM. 

We require the decay temperature to be Tdccay ^ 5 MeV, in order not to affect the success of BBN predictions, 
and Tdccay < m i) L / x f^ because the decays of r\ L happen after they decouple. Thus, the range of masses for which N\ 
could be a good WDM candidate is 

24keV (|£) < (M0 WDM < 2.4 MeV (^g^) • (32) 

Finally, in order to choose suitable sets of parameters there are a number of constraints that need to be considered. 
The null result for the process e + e _ —* Z* — > tjhVl in LEP II searches for neutralinos, imposes the bound m riH > 
120 GeV when m r)L < 80 GeV (48j |. Alternatively, in a range of parameters we will not explore, the neutral inert 
boson mass difference must be m r]H — m VL < 8 GeV [48[ for m, lL + m, lH > mz due to the LEP I measurement of the 
Z-width, which implies m nL > 40 GeV. In addition, the suitable set of parameters should also be within the allowed 
range provided by electroweak precision measurements [26l . |32j . 

There are also constraints on the A couplings in the scalar potential, Eq. [51 Vacuum stability of the scalar potential 
imposes (2(| 

Ai, 2 > 0, A 2 < 1 

A 3 ,A L - A 5 - |A 5 | > -2y/\J^, (33) 
and perturbativity of the scalar potential imposes [26| 

A 3 2 + (A x - A 5 ) 2 + A 5 2 < 12A! 2 . (34) 



8 




100 



m [GeV] 

1l 



FIG. 5: Same as in Fig. H but for Mm ggs = 500 GeV, m VH = 150 GeV and m ± = 300 GeV. 



In Figs. 2] and [5] we show regions of the m VL — Mi plane in which Ni has the right dark matter density for two 
different sets of parameters. The Higgs mass is MHi ggs = 160 GeV, m VH = 125 GeV and m^i = 130 GeV in Fig.[4]and 
the Higgs mass is A/mggs = 500 GeV, m, lH — 150 GeV and m n ± = 300 GeV in Fig. [5] The upper panels of the figures 
show the bounds on Ai, obtained from vaccum stability (cross-hatched violet regions) and pertubativity (shaded gray 
region) arguments. From top to bottom, the blue, red and green colored narrow bands in the lower panels of Figs.0] 
and [5] show the regions in the m VL - Mi plane in which flN 1 h 2 in Eq. [28] is within the 3er measured range for the dark 
matter (either CDM or WDM). The different colors of the narrow bands indicate different values of Xl, as shown in 
the panels. The unshaded background region labeled "WDM Possible" corresponds to the range in Eq. (2U where the 
Ni may constitute WDM (above it, it can only be CDM and below it, only HDM). For any particular value of Td eC ay 
between 5 MeV (which defines the upper boundary of the unshaded region) and m r)L jxf (which defines the lower 
boundary of the unshaded region) there is one value of Mi given by Eq. 1311 within the unshaded background region 
for which Ni would be WDM (Ni would be CDM for all larger values of Mi and HDM for all smaller ones). Notice 
that the lower boundary of the unshaded region depends on A^ through Xf, thus the blue, red, green colors of the 
lower regions, for which the Ni can only be WDM. Thus, within the unshaded background region Ni could be WDM 
or CDM, depending on T^ ecay . For a given set of parameters defining the model (and hence a given Td eC ay), in order 
for iVi to be allowed as HDM, its mass Mi must be, at least, a factor of ^DM^/^HDM-max^ 2 = 0.1099/0.014 ~ 8 
smaller than the value at center of the colored band defined by Eq. [28] for a given m VL . The figures show that the N± 
could be HDM even for masses as large as ~ 1 keV. 

In conclusion, we have shown that Inert-Sterile neutrinos, produced non-thermally in the early Universe, could be 
a viable WDM or CDM candidate. They are virtually non-detectable in either direct or indirect dark matter searches 
because of their extremely weak couplings to SM particles. Thus, their existence could be revealed only by discovering 
other particles of the model in collider experiments. We should keep in mind that the dark matter may consist of 
an admixture of different types of particles and that particles undetectable in dark matter searches may be part of 
it. The existence of these particles could only be inferred from collider data, supplemented by the null results from 
dark matter searches or with results from these searches which find other detectable dark matter components with 
a density smaller than required to constitute the whole of the dark matter. Unveiling the nature of the dark matter 
does necessarily require the combination of collider and direct and indirect searches. 
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